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ABSTRACT 


A fundamental, non-linear, ordinary differential equation is 
derived for the flexural vibrations of a non-uniform beam, including 
shear. 

This equation is solved with the specializing assumptions of a 
constant shear coefficient and that the moment of inertia is directly 
proportional to the cross sectional area. 

A perturbation method is used and the eigen frequencies and mode 
shapes obtained are functions of a perturbation parameter ¢€. These 


solutions are accurate up to the second order of e. 
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SYMBOLS AND UNITS 


ENGLISH LETTER SYMBOLS 


A 


A* 


Bx 


Beam cross sectional area, ie 


Lengthy combinations of other quantities (see equations (14) 
and (14a) 


Proportionality constant, I = CA, te 


Young's modulus, FL. 
Shear modulus, FL? 
Moment of inertia, te 

Shear coefficient, dimensionless 
Beam length, L 

Bending moment, FL 

kAG(y'-}), F 

Time, T 

Shearing force, F 


Distance along beam, L 


Beam deflection (positive downward), L 


GREEK LETTER SYMBOLS 


Ss 4G yee 


qx : 
, radians 
EC 
kKG@> L 
, eee 
Mass density, FT L 


a] 
kAG, F 


4 


Slope of beam centerline due to bending 


Frequency of vibration, ie 


OTHER SYMBOLS 
/ 


the 


Symbols not seen here may be found in Appendix A. 
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INTRODUCTION 


The Timoshenko beam equations have been well established [5]*. 


Neglecting rotary inertia, these equations become: 
s-ler dy] +kAG(d-¥)=0 (1) 
rAd - £[kKAG(H-Y)]FO (2) 


The left hand sides of both equations (1) and (2) are functions of 

y and 4. Because of this common dependence, these equations are 
referred to as coupled. Uncoupling of equations (1) and (2) involves 
their reduction to one equation in which either y or t, but not both, 
appear. 

It has been shown that the correction for shear is much greater 
than that for rotary inertia [10]. The effect of shear becomes 
increasingly important as the frequency of vibration is increased 
[9]. Therefore, the effects of shear are worth investigating. 

Figure (1) is a free body diagram of a beam element. Note that 
the quantity (dy - pr) in Figure (1) is the slope of the beam 


centerline due to shear deformation. 


Figure l 





*List of references appears on page 53 


The shear term in equation (1) contains the shear coefficient, k. 
Timoshenko defined the shear coefficient as average shear stress over 
maximum shear stress. This shear coefficient, then, is a function of 
the cross sectional area of the beam. Timoshenko found it to be con- 
stant for general classes of geometrical shapes. However, there have 
been other studies investigating this shear coefficient which have 
disagreed with Timoshenko's values. In 1966, G. R. Cowper [4] 
derived expressions for the shear coefficient which were functions 
of the geometry, Poisson's ratio, and the dimensions of the cross 
section in some instances. R. D. Mindlin and H. Deresiewicz [9] 
maintain that Timoshenko's shear coefficient is a function also of 
the frequency of vibration. The following table contains some 
representative values for the shear coefficient. These values 
were determined by different methods and different assumptions. 
Therefore their comparison is difficult. 


Table 1 (From Cowper [4]) 


k for rectangle k.. fier cinede 
Timoshenko 0.667 Or 0 
Mindlin O822 0.847 
Cowper v = 0 Owes 0.857 
De= ed OF S50 0.886 
v= ,5 0.870 o22700 


There seems to be a conflict of opinions in the literature on the 
shear coefficient. This paper develops a general equation with k 
as a function of beam length but all solutions treat k as a con- 


stant. This is justified by the fact that for the class of areas 
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investigated, that is,where the moment of inertia is directly 
proportional to the cross sectional area, there is general agreement 
that k is a constant. 

The remainder of this paper will be devoted to the uncoupling of 
equations (1) and (2) and the subsequent solving of the uncoupled 
equation for its eigen frequencies and associated mode shapes. This 
uncoupled equation will be general with respect to cross sectional 
aveay moment of imertia, and shear coefficient. That is, these 
quantities will be functions of the distance along the beam. 
Solutions of this type would be of use in optimization of design. 

It is often desirable to minimize deflection, bending stress or modify 
frequency of vibration in beam design. By allowing the cross 
sectional area of the beam to vary with x it may be possible to 
reduce deflection and bending stress or to modify the frequency of 
vibration through appropriate selection of an area function. In 
this thesis the area function will be defined in such a way that 
only one dimension of the cross sectional area be a function of x. 
For the sake of comparison to existing solutions, this variable 
dimension must vary in such a way as to keep the beam volume con- 
stant. That is, the beam volume must be the same if the variable 
dimension is a constant or a function of x. This stipulation will 


become evident in the derivation of the area function. 
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DEVELOPMENT OF THE GOVERNING EQUATION 


WLEh che Subst Eubions: 


S=kAG(y/-¥) (5) 
= —|_. 
ar AG (6) 
equations (1) and (2) become: 
- / 
(EIy| +s=O (7) 
fAY =e (8) 


where ‘s and ~* = 


Solving equation (5) for Jf: 


Y= y'-so ©) 
Ye ye -(sc)’ (10) 
Um y"- so)" (11) 


Equation (7) is: 


Ely"+EIy'+s=0 (12) 


Substitution of equation (10) and (11) into (12) gives: 
Ely" s'a-Qso/J+ET[y"- s’o]+ s[\-(EIe’)'|= © (3) 
Substituting equation (8) into (13): 
ET Ly" SGAY)-20 PAYPET [y- GAY] + s[-(EIeyJ20 04) 
Writing equation (14) as: 


Mir sB'=@ (14a) 


S= 7 (ts) 


Be 
where Ee (see Appendix B) 


(We 


Differentiating equation (15) with respect to x gives: 


a } ae At’ T+ 
pAY = = Bre 


pAy Ry At at REO (16) 

This form of the equation does not seem to have been previously 
reported. [It is clear that this single equation is fundamental to 
the study of vibrations of non-uniform beams, including shear. 

For the case of a prismatic beam equation (16) reduces to: 

Ue 4 2 

FI Sy, - Eira. + PAS =O (16a) 
From this point on, the analysis in this thesis will focus upon and 
be confined to the case of simply supported beams for which the end 


conditions are: 
(Wiac a (Y dt = [Puce =bW i=, =O (17) 


These conditions prescribe zero deflection and zero moment at both 
ends of the beam. Solutions to equation (l6a) are: 
, TT 
giw= Sindnx » ),= 42 (19) 
where: ub 
~~ PEre ya 
Leef dn + pA 


are the eigen frequencies of vibration for the uniform beam. The 
associated mode shapes are the corresponding gis: 

Equation (16) is completely general and it is clear that difficulty 
will attend attempts to find useful solutions which also satisfy 
equation (17). Therefore, in order to make the equation more tract- 
able, a specialization is made at this point. The area function is 


limited to widthwise variations so that the moment of inertia is 


1a} 


directly proportional to the cross sectional area and the shape of 
the cross section is Limited to those for which the shear coefficient 
is constant. 


After expansion and combination, equation (16) becomes: 


ALECy™ yi pw'y +Oew* yil+A’ BEG yi —- 
ees Yoeatysaceyi]- 4 KOsury! 
+OEC, yi A Tot eee ian -Geury | + 
ja gn Tae "4 SC pw ‘| i 
CANA’ aA As aK BE [ogee yi oF an 
ti leccyrsser a 2y'] - 5 Baas - AA’ Sat aa vo 
+4 9%) i] + AMAA Ede 
Sp Ry] - Fad a2 Tec yr -WBpw'y +8 rey 
+AVECy"+ Spo" y] cei SoA A AN 
[SS’arry] + AS Ye roury| + AX ASA 
l2e* PU zy ( = 


Equation (20) is the governing differential equation for the 


(20) 


vibrating Timoshenko beam, neglecting rotary inertia, when the 
shear coefficient is constant and the moment of inertia is directly 


proportional to the cross sectional area. 
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SOLUTION OF AREA FUNCTION 


Because of the high degree of non-linearity of equation (20), it 
does not lend itself to closed form solution. Accordingly, an 
approximate method of solution is employed. It is desirable to 
obtain a solution which permits verification for a prismatic beam. 
With these ends in mind, a perturbation expansion of the symmetric 
mode shapes is considered. 


y= Sm Tx + €SiINSTX + Ee" sin STpx on 


Physical arguments suggest that the substitution of equation (21) 
into (20) will generate a symmetric area function. Substitution of 
equation (21) into (20) yields a differential equation with the 
parameter ¢€ in it, and we will take e€ to be small. That is, a 
perturbation method of solution will be employed. This solution 
will be a perturbation from the fundamental ade of the prismatic 
beam solution. If the variable area function is appropriately 
defined, the perturbation method will provide a link to the 


prismatic beam solution. This link will be effected by setting 


the perturbation parameter equal to zero in the expression for 


the variable area. Let: 
A=A.,+€A, + €7Az (22) 
ors We +e Ud? 4 €* Wwe (23) 


In equation (22) Ay is presumed to be constant since the solution 
We are seeking is a perturbation with respect to the solution for 


a prismatic beam. 


LS 


It is expected that the solutions of equation (20) will be 
valid only for small values of e€ since they will be perturbations 
from the prismatic beam solutions. 

When the perturbation expansions (21), (22) and (23) are sub- 
stituted into equation (20), area terms of higher order than ae 
are generated. These terms may be neglected because the perturba- 
tion, in this case, considers only terms of second or lower order. 
For example, consider the third term of equation (20): 

AA's [eA +e? A.) [e Al + 72 AZ] 

: [e7A, + et AT +2 2A, A, |e Al+ Az 

The lowest order of ¢ generated here is sou therefore, all terms are 
of higher order than wl and the term hoe A' may be neglected. 

When equations (21), (22) and (23) are substituted into equation 
(20) and common powers of e€ are equated, three linear differential 


equations in A A and A, are generated. The equation for rg is 


ies 2 


simply the constant cross sectional area case. 
EAC yo + OD Aap ye! = PAWS Y= O 


For a simply supported beam: 


er _ ECE) 
°" p+ Orth 


Do is the leading frequency term in the perturbation expansion (23). 
The first order of e generates a linear differential equation 
in A, (x) 
Ai ES ye - GEC, ye"] -Al [ROpus'y, - Spud ye! 
“ECiyd’- CECE] + AL DEC y+ A, (2EC y= 
-lp Wi y+ LOpuUs yl |= -A, [EC y™- pos y, 
“PU? Yo + Qe ws yi + Selo ye] 


This is a linear differential equation with known variable coefficients. 


(24) 
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The only restriction on the total area function is that the volume 


of the beam must remain constant. Therefore: 


\A dx= Volumes A.L+€( Ads a e*( Aadx ae 


The perturbation parameter e¢ is arbitrary thus defining a class of 
areas. Therefore, to keep the beam volume constant for all ¢«, it 


is sufficient to require that: 
\ A, dx = \ Ardx=O (26) 


Assuming that A, (x) is piecewise regular and symmetric it may 
be represented by a Fourier series of the form: 
o©® 
AS > [Nesin pies +Micos TEX | 
= 
Equation (26) restricts the arguments of the Fourier expansion to 


even powers of ' 


(A, dx= 4 > [Nesin em +Micos (ix ]dx=O 
~ os [- dr Ni (cos iw-1) + fe Mi sin ar] = O 


Some, (= ©) 154, 3--- 


Cosim-1= 0, b= 0,23,4,-.. 
Therefore, for ( A, dx= © , it is sufficient that the arguments 
Q 


Of the Sine and thesecosine fumctions in the Fourier series be evens. 


Substituting: 


No= Sin Ts 
Pasi By 


into equation (11) gives: 
til if 
A, K,cosx +A, Kzsina -A, K, cose -A, Ky sine= A|K.sin2« +k, Sina (27) 


Examination of equation (27) shows that only the cosine terms 
in the Fourier expansion remain in the solution for A, (x) ; it gene 


sine terms were allowed to remain they would generate products of 


17 


cos @ and cos nQ and products of sin a and sin na. These products 

can only be represented as cosine functions. The right-hand side 

of equation (27) contains only sine terms. Therefore the coefficients 
of the sine terms in the Fourier expansion must all be zero. On the 
other hand the cosine terms in the Fourier expansion generate only 
products of sine and cosine which may be represented entirely as sine 


functions through trigonometric identities. Therefore: 


Ge 
A) = > Mn COs no 
n=O 


Substitution of this expansion into equation (27) gives: 


> M[te? K, 4 sin(ne+<) +4 Sin(na=a)t - pry kK, \4 sin(natat) 
+t sin(a- not) (eX) k. i SINCnox +o) +4 Sin (nx—)} — (28) 


Ky  tsin(nata +4 sin(- rage | = Aolks Sin 3% +k; sin | 


There are only two values of n which can satisfy this equation: 
A==2 
n is taken equal to + 2 because, 
COS At = COS (-Qe) 
Substituting n = 2 into equation (28) gives: 
M.Z,sing%+M, Za, sink = Aj|kesin Ba +Ke sina] 
Expansion of the terms K, and Z5 show them to be identically zero 
(see Appendix A). Therefore, equating coefficients of common 


arguments of the sine function gives: 


M= SA, 


A,m=M,cos ATX (29) 
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The unknown frequency is contained in the term Kg. 


MEIZ SS ES, Ka G) 
= A.ew? [I+OEY =O 
Therefore; 


Go =e (30) 


Thus, the first eigen frequency has no contribution from the first 
order e solution. 
The second order terms of ¢ in equation (20) produce a linear 
differential equation in A, (x) with known variable coefficients. 
A. K,cos« + Ab Ki sinx- A, K,cosx -A2 Ky sine = 
FISK + Fy Sin Sx+ Fa sin Su (31) 


Equation (31) is similar to equation (27). To maintain constant 


volume, it has been required that 
L 
( A,«)dx=O 
0 
The solution of A, (x) , therefore, will be of the form: 


A.) = S M, cos now (32) 
=O 


Substitution into the differential equation (31) gives 


The value of n = + 2, +4 are taken because cos x = cos (-x). 
A) = Mzcos 24+Ma,cos 4a (33) 


This is the part of the area function whose coefficient is ¢ 


Substitution of equation (33) into equation (31) gives: 


MaZasind+[M3Z,+M3Z3 (Sin 3% +MaZy sin Sa= 
B, SIN + Fasin Za + Sin sx (34) 


is identically zero (see Appendix A), therefore: 


M.Z,= eS 


Zo 
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The frequency «,* LSesconCainedsin Chemeermn © 


€.= 4 - Acwilp+er El l= 0 


a (35) 
“a = “AT + Or Ell 


To find the unknown amplitudes, coefficients of sin 5a” and sin 34 


are equated. 


ME ae 
M.= s.- MsZs 
Za 
At this point the class of areas in terms of the parameter eé is 


: Di : : 
determmened to the second™order cr ¢. Do 1S tne Contribution to 


2 
the first eigen frequency from the e€ solution. 


A =A,+€ M, COS Jo + €*{M,cos 2x +Mzcos 4x (36) 


The frequency of the first mode of vibration has also been deter- 


mined to the second order of e. 
a_ Sqaajue en 7 
Pli+GEs] 


Equations (27) and (31) are third order ordinary differential 


Uo (37) 


equations. They both have three boundary conditions imposed on 


them. For equation (27) the boundary conditions are: 


x=O |) A=M, 


) 
X=L ‘ A.=M, 
L 
(A dx=O 


The conditions on A, are™ 


X= O 4 A,= Mat Ma 


ee : A,= M,+Ms 


JA, dx=O 
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The solutions obtained for A, and Ay satiieah@y both the differential 
equation and the boundary conditions in both cases. These solutions 
are then, unique [13]. 

The class of areas obtained here are symmetric about the midpoint 


of the beam for all values of ce. For e = 0 it may be seen that the 


solution for the prismatic beam is obtained. 


Za 


SOLUTION OF EIGEN FREQUENCIES AND MODE SHAPES 


Extensive algebraic manipulations have been carried out through- 
out the following work. For the sake of clarity of presentation, many 
long algebraic expressions are represented by symbols. The expansion 
of these symbols may be found in Appendix A. 


Using the results obtained for A, and Ay » it is now possible to 


1 


obtain equation (20) in the form: 
Geer fmedyn +L Oneriyy +hine yi +eeweyeO 38) 


For the simply supported beam, the boundary conditions are: 
IS NS ee © 
ye lL )Y=O Nea oO 
The procedure in solving equation (38) will be similar to the 
procedure used in obtaining the variable area functions A, and A,. 


I 2 


That is, three separate differential equations in ea aE and 
2 will be obtained by equating like powers of e on either side 
of the equation. Solutions to these three equations will give all 
the eigen frequencies and deflection modes in terms of x and the 
perturbation parameter e. 
The notation used will be as follows: 
We > Sino Ss Ynn * © re (39) 
COs = CO CONE (40) 
The first of the double subscripts refers to the eigen frequency 


or deflection mode. The second subscript refers to the order of ¢€. 
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The zero order of € gives the constant cross section equation. 
IZ a ~ 
EC Yna +2 P20 Yno — PUn2 Yno = O (41) 
2 Z Wt a 4) _ 
(SD ac ii (Vie) v (Sind) eo Cn X= © 


The mode shapes are of the form: 


ies = Bno Sin Anx : An= a= (42) 


The amplitude, B , may be determined from initial time conditions. 
Substitution of equation (42) into equation (41) gives: 

Bro Sin Ax[EC, AP-Gewa Ah- ewe |= O 
Therefore: " 

Lj, = EC, An (43) 

eLi+sS Me 

This result, for n = 1, checks with the earlier result in Section 
ay . 

Equation (24) may now be expanded into a differential equation 
im. Yai’ By equating the coefficient of e¢ in equation (38) to zero, 


We get: 


ye + Qa ye! Qs Yar = Qy SIM NL HQe Sind COS NA 


—-Qe cos Re SIN NX 


Applying trigonometric identities: 


Qn + Qa Yai -Qs Yn > Cy SIN Nx (44) 
~41Qs4+Qe[sindn-are +4 [Qs - Qc] sin (nt2 bx 


where, 
eS a 2. 
Qy= =e eLl+e An On 
Tf Q), were non-zero, its contribution to the solution of equation 
(29) would be non-periodic [7]. Therefore, its value must be 


chosen such that it vanishes for all values of n. 


Z3 


A Fourier series is assumed for ee ee The coefficients of all the 
cosine terms in this series must be zero because the right hand side 
of equation (44) contains only sine terms. After substitution of the 
Fourier sine series into equation (29), only three sine terms have 


non-zero coefficients, 


Nat = ny Sin(n-2)oe +i Sin ne + Dri Sin (na) ol (45) 


The requirement of Q), = 0 for all n makes Cl +" forall n- 
Substitution of the assumed shape for Val into the differential 


equation (44) gives: 


BulQ pt Qa PQs [Sin(n-d)a + DrilQqt-@, 4g Q: | 
Sin(nta}on = Qy Sino +4HK),-Q, [5 n(n+r)a -4)Q.+Q[Sinbe-ala 


where, 


p=Q-ayt 


q= (n+2) = 


When n = 1, the sin(n-2)q term becomes -sin QQ. 
Therefore: 
Qu - Ss - Ze | =O 
Expansion of the quantity 4(0+@| for n = 1 shows that it is 
identically zero (see Appendix A). Therefore: 
Qy = O 
Bro PLIFOE wi =O 
Wii ZO 


2 
This confirms the earlier finding in the A, (x) solution of w, = 0. 
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To maintain a periodic solution, Q), must equal zero for all 


Gther vatues of n. Therefore: 


A 
The amplitudes Bul and Di are found by equating terms with 
the same sine function arguments. For n = lL: 


3,,= [Qy - aE Qc] 
@, pi Qape as 


ee ee 
Du= 2iQataaq- Oc 


The expansion of Day shows it to be identically equal to unity. 


This checks the earlier substitution of: 


Y=Yor EY +e? yn , Y= dsnm3x, n= 


Forwother'n: 


Bu = —z@s+G@) n> 


212 p+-Qas p2-Q; | 


Table 2 







Amplitudes Frequency 


> ae i. 
2IQ, P*t-Qz p*Qz Zia, gt-Qs 9°-Q, | 


2 e e e ° e 
ay 1 tswEne COheribDutELON OL Ene first orderve solution to the overall 
n 








; 2 
eigen frequency w.. 
a mel BE ot 
Ln = Wo tE Way + E7 Wrz 
Bl and Divi are the amplitudes of the sine functions which 


contribute to the overall equation for the mode shapes. 


Yn= Ynot © Yn + €* Yang 
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where, 


\fr = Br, sin(n-aja SED) s, sinCnrt a\~x< 


The second order ¢ equation in the unknown y 9 is; 
n 
oe 4 


+S3 sin(nt4)a + Sy sindn-4)a + ds sin na (46) 
The assumed shape of Yn2 is therefore: 


Nna= Ina Sin(n+ad)a+ Dna SIN(N-A)x + Ena Sin(n+4)og 


+ Fra SINCA“4H) XK + Gar SIN RA 


The term O. is. 


ds= Bao PLI+O An | Wan + Se 


As was the case with the first order ¢ solution, this term may not 
be allowed to persist if the solution is to be periodic. This 


condition defines the frequency and demands that G5 = 0. 


For n = l, the 6, sin(n-2)q term becomes - 55 sina. 
Therefore: 
ut = d2-d¢ 
_——- a 
i Rie eC +O()*| As 


Similarly for n = 2, the by, sin(n-4)qa term becomes - Oo, sin 2a 


and combines with the oO. sin 2a term. Therefore: 


oe Sy - Se 
“a2 Bao pLI+4OE)2] Az 


The amplitudes B E and i are found, as in the first 


iy 2 


order e¢ solution, by equating coefficients of the sine functions 
with like arguments. Care must be taken, however, noting combina- 
tions for certain values of n. Substitution of the assumed form 


of y , into the differential equation (46) gives: 
n 
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Bna lq! -Qaq- Qs] sin(nt2)« + Dna [Q, BtaQy p*-Qs | 
SiN(n-alx + EnalQqr*- Qagq* Qa sin nada + 

FaalQ, pt Qap't Qa] Sin(n-4) x =]S, sin(hga}d ¢SzSin(tea) oe 
+Sssin(nt4o+ Sy sintn-qa+ Ss sin na] ze 


where, 

p‘=hn-4)E 

qe) ZT 
For n = Ll , 
fa,GZ)+ @, G)+ Qs ][Bia- Fa] =6.- du) ae 
The expansion of 5, - 5, shows it to be identically zero for n = 1 
eee AppendixwANh elthe expassadion QG)'-o.)- Qs | ieunon@a@eo. 
Therefore [Bi5 ~ Foo], which is the coefficient of sin 3a for n= 1, 
is zero. 

DialQ,CE}"- QacEl*- Qs] = Ss - $2) £2 
From the earlier calculation of ie we know that [5. - 55 | = 0. 


12? 


Therefore: 


Diaz © 
E.2[Q,2)'- Q,(&2)*- Q, | = $3 /p2 
Expansion of “Yagi Q97-Q3]A; shows it to be identically 


equal to one for n = 1 (see Appendix A). 
E.2.= | 
This checks the earlier substitution, the the Ay solution, of: 


y= Yor EY, +E7Y, » Ya 4sin OTs. n=l 
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; 2 ; 
Expansion of the expression for 5 also shows it to be consistent 


for the value w, found in the Ay solution. 


For n = 2 the sin(n-2)q term vanishes. Similarly, for n = 4 the 


sin(n-4)qa term vanishes. Therefore: 


DE = O 
Fuaz= O 


Forvall “ether in. 


= F-yi SIeIK n> | 
Ona Qi qt-Q2g*- Qs As 


5 8 ee, 
Dna= [Q.p4- Qa p*-@el AS 


=e a re ey 
na Q.qr4 = Q:iq/*- QalA. ¥ 


— yy ny, | 
»  LQuP*- Qap- Qaj Ao 


W) 9 1s the coneribution of the =seeond order © solution tovtnhe overateL 
n 


: 2 
eigen frequency W . 


aL y a Fi 
OO (tO, Cae 


B 2? D2? E42? and P 2? are the amplitudes of the sine functions 


which contribute to the overall equation for the mode shapes. 


Yn= Yno a ENai €* Vina 


where, 
\na= Bra SINGH AK + Daa Sinead) Ena Sin(nt4 Jo + Fra SINCN-Y) of 


The mode shapes for the class of areas derived in Section IV are: 


n= BroSin atx +€/ Basin PX + Dai Singx] + 


"(Bra Sin qx + Dna SiN px +En2 Sing’ +Fry sin p’ | (47) 
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PO - fh) - (a)'O]) [O-LxO-(WOT] [FO -GyO- (AO 








a7 [ (HOt) ]d EG or 
Barz 





Ce lysed 


The eigen frequencies are: 
= Wy 
i 2 Brno EG: An Rye 


= (48) 
é BroeLit©O AR 


where, 


Xe (Sa-Scne) - 

Xt (Sy -Solnna ka 

oe = ~ Se x2 5m >or 

oO is an amplitude which depends on the initial time conditions 
of the problem. 

By setting ¢ = 0 in equations (47) and (48), the solution is 


reduced to the constant cross section solution. This is expected 


because when ¢e« = 0: 
A=A,= constant 


By setting n = 1 in equations (47) and (48), the values of 


and 0, are 


a 
= Sin TX in Stk + Sin Sw 
Wi= S TK + ES = = te 
aes Zz 
a EC (2) +€ x , 
UO), = a oe 
eLi+Oe)?] 
These resultspare tdentical with the results andesubstitutions from 
Section IV, the solution for the area function. This checks the 


great amount of algebra which was performed in obtaining formulas 


for the eigen frequencies and mode shapes. 
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CONCLUSIONS 


A fundamental non-linear ordinary differential equation, equation 
(16), has been derived for the flexural vibrations of a non-uniform 
beam, including shear. 

Equation (20) is the particular form taken by equation (16) under 
the specializing assumptions that moment of inertia is directly pro- 
portional to cross sectional area and that the shear coefficient is 
constant. 

Tables 2 and 3 contain approximate solutions of equation (20) for 
all eigen frequencies and mode shapes in terms of a perturbation 
parameter ¢€. The perturbation solutions are accurate up to the 


second order of e. 


ot 


EXTENS LONS 


The preceding results could be applied to beam optimization design. 
Using the formulas of the preceding sections, the deflection may be 
obtained as: 

y= FCx,t ,€) 
It may be desirable to minimize the deflection y. For example, 
minimization of the amplitude of vibration could cut down on noise 
level. Bending stress is a function of the deflection and this 
quantity could also be minimized. This minimization process could 
be approached through iteration on the perturbation parameter e. 
A value of e would be chosen and y(x,t,e) would be evaluated at 
increments over the beam length over appropriate time intervals. 
The maximum value obtained from this process would be saved. This 
process would then be repeated for other values of ¢€. The minimum 
value of the maximums noted for each iteration would give the optimum 
value for e. 

It must be noted that the equations developed are valid only for 
small values of e. 

Another possibility for further study would be the application 
of this method to the forced vibration case and subsequent optimization 
of design for deflection, stress, or frequency. 

An experimental investigation to compare with theoretical results 
would be useful. Problems would probably arise in the manufacture of 
of the beam itself and the development of devices sensitive and accurate 


enough to detect variations from the solution for a prismatic beam. 
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APPENDIX A 


EXPANSION OF SYMBOLS 
This Appendix contains the expansion of the symbols used in the 
main body of the thesis and those symbols used in Appendix D. Symbols 
equal to zero are shown. 


The expansion of the symbols used in the area solutions follows. 


kK, = O7eu2 @) -SEC EY 

K,- 3@pw2 +S eUsE)+ EC (@y-OEC(E)* 
k= -2EC (EE) 

k= 2[ECKE*- ews - Se WHEL EO 

k= EC CZ)t- pws - Sew. (=)* 

Ke -GewFGZP-ewR=O  , (W,*=0) 

Ky= - PW - Sew (E)* 

Kez ~Pwr - pw? @t)*= © 

Kaz (=C,(S2)* - eros - @e w?F Er)* 

Kz 2KQ=O 

Knz -2@"ewe- SEC (E)* 

Kinz 77 pwe(t) -GSEC, CE)” 

k= @7ewa (E)*- OS EC(E)* + 7 ew? 

kw= - BECKY - e*e TOE: 

Kez -B@ pio? - O*pwiEF =O 

k yz QEC,(RE)* - BS pu? -S% ews (a)*- EC, AZ)” 
k= O?p we (AZ) -S EC, (AEP 

Kig= O*p wPGE)=O 


B32 


Z,= +k (E) -2ka(€)+k3E1 
(ee 4K, or +2 Ky (E)“+ K2(2) = © 
Z = BAKE) + 8K (E)7+2K, (EF) 
Luz 3AK, (EE) - 8 KE + 2K.) 
= A. Ky 
G.= AoKe =O 
Cs= A, Ka 
Giz NiAo[Kio-@E)* Kis}. N= SE 
Os= NV, A.[2Ks - (ar)* Kie | 
©. = NiAo [ 54CE) Ka +(3x)° Kia | 
Or = NAL AE) Ki2 - QE)K, -@E) Ka 
Ce = NPA (2=)*K, +(4=)° ay 
O.= Nr Aol (@)" Kis +(ax)‘ Kia | 
Brox Ni AoL( amy Kiel|=O 
€ = t[QG,-Gu+@s + Gc + Gv +23] 
§.= t12GtGu+Gy -Cx% + Cy, | 
$= & R@s+@s - Ge + Cry +Grg- Cx] 
¢,= -4[Gu-@s -@c- Ge - Ga] 
Mi= & Ac 
M. = a7 MzZ2 
a a, 


Ms = ¥Zq 


The following symbols were used in the solution for the eigen 


frequencies and mode shapes. 


QQ) = aC, 

Qa = Qewrs 

Qs = Po 

Qy=Brole +Se rrluri =O tm =o , d= ne 
Ws = Pho N,(2=) ater ACE) - SEC, iy @r)*_ AEC, 2 | 
Qe= Bro any [Be punt -~6€C, Ang An +EC Ne] 
GQ)4 = Brno LEC, An -PUUns -Se Ore An lO 

Qe =Bno [2 EC, val 

Qa = Bro LBO@pu?, -SEC M4 S7evrte Ke +EC, AY] 
Wie = Bro [Sec, rs = S*e ele An 

Qu = Bro Le+ Se An aca 

Oi2= Bri Le+ 2e p? | Ce, =D eae p= (n-ay = 
@Qrz= Dai Le+ Se q’. Wi=O ge (n+ye 
4 = Eni LEC, Pp - CUur's * Deore ofl 

Qis= Da LEC. q' = eho 7 Qe CWro qr 

Vie = Bno Le+ Fe re I vom = @ 

Qit = ri 2 EC, p* | 

Cy = Dai L2 BC, o? | 

ig = Gy 

COEF = Bho 2 CS ewite +SEC, we | 

Qai = Ibn met S*ew ro Re S QE, 121 

COaa= ene SEC, CS ¥ Seu Q i + /Ge 17%, | 


3D 


GQaat= Bno [OEC, Des +O*euwr | 

Qay = Bai lOEC, a - ZG evure -6"e Core p* - EC, p24 
Qas = Dini BEC. g4- 2Gew2, - SP wie q*- EC, qi] 
Qag = (Sno [2Seumn + S'¢ Uh AR | = @, 

Qan ad Cs 

Qaz= Bni [- OEC, p*+ Oe Wrs pl 

Qas= Dn L- CEC, a ore Ure ql 

Q3o= ee [o*e Cnt = © 

Qa ee Gio 

O, =- Ao Ma [(AE)* Qe] 

b= - Ao Ms [(HE)* Qa 

b.= - Ao M3 (A=) Qs H2Z)"Qiel 

hh, = by Ma [(BQe +(4Z)FQro | 

bo= AsQu 

de> Ao Qa 


p = As Qi - > 

be= AZ N L2Qu -FEYQau | 

daz AtN: L2Qis - @E)"Qas] 

Pio ASN, LAE)" Qas - @Z) Qn! 

biz Aa NLGEY Qig - FE) Qe] 

Dia= As 2((2E) Qa, - (2) Ve -EY Qa | 
b= A: 210 Et. +(4E)” Qar 

= As NLC am)t Oaz - (22) Qaa 


&=+Ldit+ds+ b- | 

da= +14, - pz +e] 

bs =k [ha- du —de +n + OHA + $34 - Sy] 
ay [boa cu ~ bg = Do as Pix + afr “a wy, | 
ds a; [2 ds -Dg -deq +Dio = Di +12 +diy | 

Se= 4[- Py - dq + Dro err + hist dig | 
p'=(n-4) = 

q' = (n+4) 


APPENDIX B 


THE CASE WHEN B* IS ZERO 


On page 10 an assumption that B* was not equal to zero was made. 


Suppose B* did vanish: 
4 4 
BY=|-(EIa‘J=O 
/ 
[EIo’| =| 
Integrating with respect to x: 
dg . 
EI Se= x+a 


Lf tf = Cconstane anda [s— CjA 


en Maem 
Oe ,U = KG a 


A x 
ea dA - 
rare ( A (+ a)dx 
~ EC _ 
ee hn A : Ze +ax+ b 


A= ex p[- Se (<¥+ OX + by 


This is not a symmetric function about the beam midpoint. Since we 


are restricting consideration to situations having symmetry about 


x = L/2, this case is excluded. 
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APPENDIX C 


NUMERICAL EXAMPLE 






‘Hebi 
L= 100mm—_—_—>] k—béer-o| 
A=ActéAit+e' Ay 
y 
Steel Beam = = 20x08 os) A.= tint 
Ge1axlo°ps: Brot Lin. ny! 


k=.67 IT=GA C= hAja= Kain’ 
p = 490 bt fe 


A computer program was prepared that computes the first ten 


eigen frequencies and mode shapes for ¢ = .05. The program also 


computes and plots the variable area function eA, ts eA, fOr 
e = .01, .02, .03, .04 and .05. The first ten eigen functions for 


¢é = .05 are also plotted. 


The physical significance of cA, ts eA, in this case would be: 


A= boh= Ag+ €A.i+ e7Az 
~ A A, 2A 
nae “ee tS 


Therefore, the width would be a function of x while the height, 


h, remains constant. 
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The first ten eigen frequencies are; 


10 


2 
Ln 6 
4 


~332x10 


5 


~530x10 


-268x10° 


6 


-845x10 


7 


.206x10 


7 


4 25rG 


7 


~785x10 


8 


33°10 


8 


eoi2x10 


8 


,oee x10 


40 


.488x10 


.941x10 


.124x10 


pa 
CAN 


6 


-.261x10 


7 


- ,310x10 


909x107 


8 


. L90x10 


8 


SEPM, 


8 


8 


.693x10 


8 


9 


9 


. 5 éx 10 


UM= 


5136 "eyeles/seam 


al 
53902 


944.9 


1,462. 
2 09cm 
2,832. 
3,683. 
4642. 


5,708. 


a ele 


X Cin.) 





(24) *V 2D ey 3 





Ga.) 


(ev) Wy.3 + ® 


pepe: 60 gO 1009 


‘a | 0 ts 


(qv!) iG) 22+ Vv? 
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F, rst mode sha pe 


oO 20 40 60 gO 100 


ro) ao Ae) gO 100 
X (n.) 


ied mode shape 


——— 


SSE ie 


ot 


0 


Fourth mode shape 


Bh 
yard o 
“< 20 40 GO 80 too 
% x (IA) 
Picea mode shape 
a 
yirlo 
a 
ao 49 60 80 too 
x Gn.) 
Sixth mode shape 
ol. 
YiMo 
-a 
vo) ao 40 co gO 100 
x Cn.) 
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LOOP FOR CCMPUTING THE VARIABLE AREA 


Lib 

& 

= 

+ 

ae 

ae 

= 

= PM! 

CGC 2 
Lie x 
ee Om I 
OOet om 


ws HCl! 
mre Il + 
OO wom Ye 
t% % lt 
e@ eQwwf- 
NP NN 
lore O 
pm > CF) KK 


. 6) 


Ww 
ee 
N 
saad 
NJ 
ee 
= N 
~~, & 
~ ~ 
e fe 
r=) ca 
om) o 
ad 
od @ 
= , = 
ae ‘a CJ 
= i . ” 
~_ WY Cc UJ 
z a. e = 
lu NI ad 
+ ~ <I 
> a us > 
~ CO =) a 
Y u. - Ww 
©) — Le 
=) — - ~_ 
C= N toma Lh) 
% CJ oe e& 
_—— <I wal O 
9, 8 i — 6+ Ww =< 
wor $f Ce o™N (2a «<I 
(ZN GQ +* <a 
LL $e a 86+ =| Ta) 
+ # - WY e wi 
mmx Wo Ww - O 
trom food eer me! Bc oie z 
wer wow C1 = - UL 
ss FQ. - MN = 
Ht + —> Arter < Ww 
ee w tN HH - cr 
— ww od —- UF + ss LL. 
Mw © aq ANH = - 
—~ COre Ym tb! HN a =< 
€2 WW ~ HP 3 =e r— Us 
Yi eft OOF weed FMD ee Gs) 
ad ~— £HeU— FOOD MHA — 
Sa ames ODO) © oem UU or iW 
RON ANMeS AOMM CC Illa~ F 
SM ere MwA BUI & ee WN tela WwW oe 
el Zi || DOOO~w~w Ht ZaDe DP © 
am FC UW oe ee we we ZO 2 Wh 
Cesta h Wm teary me a) 
tee oom om om et ff pe pe Se et SR 
WY aA we Zee YX OY eZ 2 O 
Ca Owewew ee OY YKYOONCACAOA OO OO 
CwosLidtsosESsuUuUuUWoOr CVO 
out 
G - Ovo as ~~ 
— One — 


gsleLlAST 


Me Me ee Re ee he fe ae ae ce ae ie a ae ae te ee a ie a a ie eK ee a ee kK SJ ) 


€ 
C 
C 


47 


“4 ++ —_ ~ 
7 Ea | «<I = sa 
PA @e ye = a 
Ws NA! ui m4 
H 86+ + + yw 
—- KF € —~O 
~F —- N CO # 
= <y * H% © 
] <q <{ ~ fo aa | 
- tH MM ] 
qe el — NIA 
moO w+ <I * & 
+ ww e % % % -~ 
-—- +t HF x x= % 
~ OF — a # 
aI#ty & >< ae 
* LL tf) + & = 
us - <f (a -—-— fea) 
He Ht + = YY a 
Bo bom  _: www % 
me Nx NI ~~ som} =~ 
} +— ¥ — Ft YY 
-~- -—§O % IX OO —_ 
x~ «He xr Co © 
w OD —- = + & % 
<q Oa. * “Mo det. oO 
H+ HH = & ie xox 
=~ Oe Tt Oo ie | ~— + 
x~ awn + FF 2) a 
wv «6H UH *+ OO FS Or 
C © + ~ fF FH — FI 
H+ NT ¥YyY +t FF Ht om 
or e+ I N ~yY OO twm~ 
+ So ** + DM—- © *#O 
 £€+< us + OK KH es 
NI t~— + a. &*#kK = yO 
+ —O ~~. Fe Wt N wt 
~~ CH - ( willl + & | 
Lt ~-O 4 oo t= + OS 
—- dae -~ © oye = HH 
ww = & ~ *€# Noe MW THe 
od ot *%Xe -3YFf =r FF Cw 
—=— Ne -O w~ 2 CC Xe 
om ie MOO SF Un YY +2 
—- W eww -—F+ FF —-O * OF 
x te TV Ht ar WwW Uw TF ru 
~~ ot He FF ¢€C Fe ww 
COMO etm Ct LT — I + + + 
mt Vtomww 0 Ge Fe Oe COC om 
“Oetrsce ty FO ~~ KIN XY YY 
AO Oh eis Line HC Ce Te ww 


ReKRHAKRKRHEMR MH OC ow DD ot OF SZ 

YTOOCHO eme 20M ODI K-OONCHG 

Cw F&F COWOVNYOVNGD tT WU UOT WH 

eH UM HE UNH K PR ome ee I meme 

et AMINO EE ODE Hem ee tt dHs 

AaCOCEWC# OOF OCOCHKHOCOCATeEAeCoe 
a AN = 


men 
+ % 
+ 
wm em 
y 
co Gee 
cx 
am >< 
+ + 
LU iky 
Hit 

| 
NN 
+ 3 


~~ ee 


X(K)**4-R9%*0O(K)-TH*RO*O (K)#XX(K HE) 


28 


x 


+7 .*TH¥RO*O(K) ) 
(K )-TH*¥*2*RO#O(K)*BBCK )*#*2-Ex 


*3) 
K) *( THE XK (K ) €* 4-3, *THERO*O(K P-THE*®2*RO*EO(K) KXX(K BR Q-E * 


Ox — 
Mo» O 
% HO 
us 
% # O 
me eae ail = 4 
—}— 
In 


ww? wee Th 


oe *TH*¥RO*CO(K P+ THE*2*€ROXOC(K )P¥A(K #2) 
H 
H 


<czow~w 


HKR—-WWHE CO 


a sme 
ot e 
N N 
—_ ~~ = 
Cc om « 
% rn 
foe) owt a ome 
% — 
# awry 
= | aa 
“th © ee wer 
N Na oa. 
* “I+ 
Ca _ Gem 
-~ MAIS 
ot eat", 4 
am ama 
Wm +44 
iN- eu om 
on eee CP Na = 
OC “N+ = 
“tt —~ eu 
we NI NN 
Ot an | 
ADA Ht ang 
TUIVAI MH HN + I 
NN St NN 
wor we CI HOH “0 Oo 
CORGKRE™ “1 + 
KRRHOEOO KI —— 
NANG Fe tte ac 


HP He Zeng 


Ke tHE O ADE | 


Hf mw | + eww j~— 
ome em em NID mm QO 
ter LUAON HH 9% OOO | 
were ce” | Ale © Ou om 
a CU www GONNA OO 
UW+s OCOK+4¢11— 
1’ KN ENC eon 
mmf w AY H OOS S I 
TNAIOHK Eww wm 
mame ¢ Ht Hh OOOO 
——OaMhee KR + I I + 
COKKEREREN OO OS 
HHH HANNE ASN 
om on H pe HH SE wr we ee 
amet FF —-AA0ag 


cae? eee? GS) ee ee ee ee we ee 


(LotTHeA(K )**2) ) 


+TH*T **2)) 


ae 
F WtUGM 
DOQ* 
(__ 
$m 
et & 
rol ee we 
KR Wty 
ooa~ 
we we <I 
% * 
Onno 
Dm nm X 
—S~ x Ht 
~~ wee wer TO 
-“~COF 
000+ 
—-OUOUW 
=) a 


on 
i 
So 
Ha a 
ow NN 
At 
+% + 
a EE com 
My Me 
elo ww 
~~ 86 
ax 
++ + 
im ge a 
ar <x 
Ale —_ 
~-O O 
Aart + 
(Cc G 
mY 
CH 
aI 
—eF 
Q | j 
> a 
x~ w 
ed oo 
Co we 
-~%F% ** 
OO Oo 
wr 
a } ] 
ist + 
—~t 
Ot +£ 
ow ea lie 
Oreo wc 
| — e— e 
iy >< Ha Il 
we FE om FE om 
Owx wy 
Fe wer ee ee ee 
ee SE 
N~ OO 
we gol CJ 


ond eer er eee OD) EEE er ee me ee ee Hee Qe HOMO WLU ee ree ZINN OO ANG SHIN J mt et 


*Z eZZrOOOOZNZANVUGZZOSGww SrssteKwWATowW © © © © | 


eo eOm el ell e 


BOOMOVTAANMNMANM#+ OX DM 1 MOO 1 VOWZ—wH oo of WWI HKHNNNnnhNnared—- i oeVaag 
WH HUH PP ne ee nn ta bb PrN han nem ew mem a= j) Ww i Wu 
On ee ee Om mie {ltl fl ty fl ill [mom mem HNO ST OM 0 Om I] ilo Om fC 


SCNOPM OO CAUSES YIN Y OF ODO See mm MM MeO ANOS we we ww ww Y YW HOmw yy YY & 
ee td ot td EOI OUI QUO ON me Qe NANAIMO BAO SOP OD 46 test tet I) J I IO ee I) eee a 


sCOCCSCSCCECGOOSGxXKOSSGSSOE0G00000 000000000 COO msG000mCre 
ot 


48 


+ET(KD*DSIN( 


NIA wx 

+ .- 9 =——+ Pot 

+ $e —— ee aaa 

YY “t+ = 

ese? Shae Re LL 

nN > xx UW 

+ tt «KO 

YX om, ue 86 

of ~~ © mem OC. 

+ oo YI WY 

eS ~~ OQ ul 

a Ae & Ce Ww % # + 

ae x AX eM e=t 

| - OW Ze, 

~~ LL © COW 

“MY ~ - ~ ae +430 

—— x ~ ~ mM om om H 

oC — x ~~ eNO om om em Q. 

% © — — ee ~Ntl ¥y~y~nm 

oO Oo O eel eo my eOlLL YY) ad fed wee LL 

a ala G OC — - uy ON, ef wwe +> 

i | o ~~ wf an] a ~Ne2? OMY xzst>— 

~ a 2 -~ ~ _ 1 ~ ~ SS +77 OG % tum 

% % x (vr) - ar x x Ze wu LL mom OM 

5% pe on a be — ~ Ie Oo WW bt ee 

= Ce oo i = a — MODADNY wee wee C/) eee 
Ow MOIS S " = = © a) N ob Wl Ox OS 
ew= @€ 6 © ~ = ~ -_ ~ e eDnTr ca oH 
CHIYKCOC ~ - ~ = ~ ~ AY DD Se ees oe ee 
WS WH eG ~- — Lis x ~Yrnmnmm TT - Z2exru 
MUS ww Cc _ a. = Zz a te oud J Wb <f 
ae eee ee ee ae ee? ©] ed Lo @ Gee ee Oe 200. QO OU 
hm NE — rere LPWeark HK SES o 
Common lL LL LL — = a coal = m~e @€ 0 eOQ0QA Od Caer ne Ct 

‘pare g aie oem ge (\| ate ate ame (0) am oom ame CT) arm orm om (TC) are oem aoe (PS) am om ome 9d 11 AVL =X ¥ = 
MO OON OH LSMONOMOMOMOFORONODIHe mire eee ee eee YO + 
N od BRINTON AI AINA AIM ANAM NN NMAANAA AINE OCOO ont a hem et ffm 
ellJUI © ee oF RH HHH HHEH HEH SH HH SHS KR HO OMIN SA WE CO CO weer om 
BGOOVWVWPWODOODODODOODOODO DOOD OOD DDDDD000 Ow www www NY | HELO 
ML YP EE ORI Cy ee ee ce cee ree tee te tee ee ee ee tee ee ee ee ee ee ee ee ee eee ee ee ee fe ff bf be et om om it bet bet 
Comm 6 @ OLLI UU UU LU WU I LW Lt a LL tl lt ie eich at ei cl ee RIN mee 


¥ MYM SEER EEE EEE EEE EERE EE EEEE EE EES SEE TES OM YEE 


See ee ee ee ee ee hawt boned Peed rund Dent Prost Pamed Peed Pad Dad Dont rand Deed frm Dead fed band Paed Pat Pat Pett ete XY OK OKO YO Wie s 


Wer UU UCSC YELLS KCKOCOOOOCCOCMOWWNWCO 
Wea SESSSABELAZEZSSESAZSISISFSSSSSUUULUUUuUUuUbsCGOONGUG 
_ 


49 


NOAM QRMNOS 
NMS NUAIAICUAIAL 
AUT eo 


glyLlA 


y//) 


PI TITLE 2009 20 eD Caley 2y2y2_5 


P22 2 tS ESS ELSES SESS SES SESS SSS ESS SELES SSS tS SES! 


wr UL et tO YO Lom 

OWOCLOXCLKOOTO 

SCUSBLSFEUUOMO 
eo 


woe free 2 bet CO ee eet et ee OF OF) CLL 
Cx<FtOCYvOre ervoorZza 
OQOOVUBUSsSSFSUUMW 


re 
~~ 
~~ 
o 
VV) 
ui 
a) 
~ eS, 
~ a 
~ ee 
& WY = 
N\ Ww ie] 
= — ti 
& WO ie 
> amd anal < tL 
ont ol uJ oe) 
- © ie 8) a “x 
~ _ <I | ut 
fad * ess Ww m(D 
Us =i * on Cre 
mM iT Se i end LL 
= «z ad al 
ae a © =z = el 
z e o Ww —_ 1 Lu 
~ Va) | a) _ are 
Ww zs x 4 — © al 
Qa 2. o eg) © oh 
<I ww —Y) Cm aw” aoe 
— x<- & = es 5 amd ~— O 7 
a) WY Wnt We z Nm b— oa Ww 
aa Lg. 3 Li —l + = > ag 0 a9) 
—_ Ww Onto - + < © N 
= O — -F o& a a.© wo LL om WY 
UL x~O fh % = a QN~ =r uJ 
om Ut = ™ eH = UES —_— amb @ a 
“OO mm PLN O eee CT em Nwet HN 
otf} OMe One w~ O w-O O-0000e = WY 
we CUO NTO F&F bm ONNONO OO — 
tL) & & ~ Lf) UW QO mew OU ON em eee 1 —_ 
flee 0O0M00ONRN~Ke DPD OO OfrIo-0OO000~~— fe a 
one eee ewe heb OO oS UO Wl emo ew ewww © 
OY) ot LE LL ELI LL eT — Co mn LT LLL <<. - 
mp bb ERIE SSO G eee Eh ee TE LY 
oud 


® 
oO oe wO ad we oO Next G 
oO wwmN Gs w& Cow 
Jt hf Oo OG r) 
— 
WOW 


50 


APPENDIX D 


INTERMEDIATE FORMS OF EQUATIONS 

Referring to page 25 the differential equation for A,(x) is: 
A.{Az Kcost +3 Kasina -AS Kcose |= As sEC, yz - pur Wi 
-PWs\2- Pur" y, +9 PW Ye + Se we ya + Oploiy"e +2AAS 
ECiyin= PW? Yo = PUY) + Sew? yo" + Spud y,"& +AoAL 
QE $ + AALS 2EC ys" 4 Al 3-287 ww? y, + OE Cy.” e+ 
AAR IS PE y, HOO EC YS + AIMS- Se wi yi" GEC VE 
+7 Op we \ok - AVA" FOEC y," - 72d Yo + AoA, "J8EC yr 
“BO pus y, - OPW? y+ Bend y+ O% ou? yo + EG y,"€ 
+ A Al S$ OEC YF - 39 PWF Yo Se We St ECE = AoA” 
98K y+ Qeur y+ Cewi\ ¢ -A, Ai SOEC: yo" + Ou Me g 
Substitution of A(x), yo, y,, and y,, and subsequent reduction 
gives: 
ALA’ K,cosx + ASK, SiNd-AS Kacosa [= GB, sind + @,Sin 3m 
+@asin So + Gy Sinx COS Ax +@s Sin BX COS Qe - 
Ge Sin AX COS Sal + G7 COS& COS ade SINAR+ 


Ge sind Seow 0 Gq Sino Sihaala 2 Bip OSX SIN Ao 


After further reduction, using trigonometric identities, equation 


(31) is produced. 


of 


The basic differential equation for Ya? is similar to that ton 
A, (x). The only difference is that Yn? is treated as the unknown. 
After substitution of Ay> Ay» Yo? and Ya? the equation takes the 


form: 


Oy Qe +Qa2 Vina ~ Q3Yna |= p, COS xe SIN No tb2 COSY SIN Nek 
+z SIN 2x COs NX - hy SINYx COS Ne +5 Sinne +heSin px 
+64 Sin gx ~ dy COOSaed Sin px - dq COS 2x SIN qx 

+Hio SINK COSPX +h SIN AW COS]X + hz COS det SiNAal COS Nek 
+3 COS*ARK SINNA +yy Sin7AKSIN NK 


Application of trigonometric identities to this equation produces 


equation (46). 
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